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We reconsider constraints on Brans-Dicke theories arising from the requirement of successful 
Big Bang Nucleosynthesis. Such constraints typically arise by imposing that the universe be 
radiation-dominated at early times, and therefore restricting the contribution that a Brans-Dicke 
scalar could make to the energy budget of the universe. However, in this paper we show how the 
dynamics of the Brans-Dicke scalar itself can mimic a radiation-dominated kinematics, thereby 
allowing successful nucleosynthesis with a sizable contribution to the total cosmic energy density. 
In other words Newton's constant may dynamically acquire values quite different from that today, 
even though the evolution mimics radiation domination. This possibility significantly relaxes the 
existing bounds on Brans-Dicke fields, and opens the door to new possibilities for early universe 
cosmology. The necessary fine tunings required by such an arrangement are identified and discussed. 
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I. INTRODUCTION 



In the last few years, there have been a large number of different approaches to the dark energy enigma. As 
is well-known, the anisotropies of the Cosmic Microwave Background (CMB) radiation [l[ and data from type 
la Supernovae @, [H, 0], are well-fit by a cosmological constant, albeit a fine-tuned one, of which we have no 
satisfactory theoretical understanding. Another possibility is that the cosmological constant is precisely zero 
(or at most subdominant) and that a dynamical component is driving cosmic acceleration. This component 
may be a manifestation of a large-scale modification of General Relativity 0, 0, 0, H, & H3, HH, 03 > or a new 
dynamical field, such MS ir.l scalar field 
non-minimally coupled to gravity 
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It is quite natural to consider scalar fields that are 

m m m m us us us m m, is m, p . such 

fields arise both in the context of string theory, and in more general theories with extra spatial dimensions. 
These fields are typically of the Brans-Dicke (BD) type, and have been considered in a large variety of other 
cosmological contexts, such as inflation and baryogenesis 

However, there exist tight constraints on BD fields. Some of these constraints are on solar system scales, 
constraining the BD parameter to satisfy u> > 40000. Other constraints arise from cosmology itself. For example, 
in order not to spoil the success of the Big Bang Nucleosynthesis (BBN) predictions, the scalar field component 
should be subdominant to radiation and the scale factor must evolve nearly as i 1//2 . This constraint is quite well 
accepted, and a detailed analysis of the constraints on other kinds of evolutions, such as dark matter dominated 
one, can be found in [48| . 

In this paper we revisit this particular issue, considering whether a BD scalar field might play a significant 
role in cosmic evolution during BBN, while maintaining an acceptable evolution of the scale factor. We shall 
only require the BD field to drive the kinematics of the universe, making it expand as if radiation dominates. It 
should be emphasized that in our model the expansion rate (Hubble function H ) is driven by the BD field, while 
the scattering rates are determined by cross sections and abundances of the standard fields of the plasma (e ± , 
neutrinos, photons, baryons). One consequence is that the relation between temperature and time becomes an 
adjustable parameter in this model. Provided the BD field can reproduce a standard kinematics H(a), one is 
free e.g. to normalize the abundances of all the fields to the value they assume at the same H (a) in the standard 
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model. Since they do not interact with the BD field but gravitationally, the particles in the plasma cannot 
distinguish which component is driving the evolution of the universe, and they will interact among themselves 
in the same way as if the BD field is not present. As a consequence, the production rate for the light elements 
is the same and the freeze-out temperatures for the reactions are unchanged. 

The big difference between this model and other BD models during BBN (see e.g. [4{J), is that now the 
BD field, i.e Newton's constant, can be largely different from today's value at the time of BBN. This is quite 
different from a standard approach in which the BD field is allowed to be just slightly different from today's 
value. However, just after BBN, some non-trivial dynamics is required to bring the BD fields to values consistent 
with a standard evolution for the universe. This is sufficient to ensure that the temperature of the universe 
has the usual behavior from the end of BBN up to now. We shall see that such a scenario is possible, while 
remaining consistent with solar system constraints. However, as one might expect, this interesting possibility 
can only arise under finely tuned conditions. We describe the tuning required, both on the initial conditions of 
the scalar and on its associated potential. 



II. SCALAR TENSOR THEORIES 



We consider the BD Lagrangian density with a potential W((p) as 
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where R is the Ricci scalar, ip is a real scalar field with units of [mass] 2 , to is the BD parameter and C m denotes 
the matter Lagrangian density. Our signature is — +++. Typically, as long as the potential is too weak to 
confine the scalar field, precision measurements of the timing of signals from the Cassini mission yield the 
bound uj > 40000 |50j . Although it has been suggested that this bound may not hold on cosmological scales 
[5ll ]. we will take a conservative view and demand that it be satisfied. 
Varying the action with respect to the metric tensor gives 
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and combining this result with the equation obtained by varying the action with respect to tp one obtains 
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where T = T/ and W >lfi denotes dW/dip. 

To study the cosmological aspects of this theory in the framework of a homogeneous and isotropic universe, 
we use the Friedmann, Robertson- Walker (FRW) ansatz for the metric 



ds 2 



a 2 (r])(-dr) 2 +dx 2 



(4) 



where 0(77) is the scale factor and 77 is conformal time, related to the usual cosmic time through a(rf)drj = dt. 
We have assumed a zero curvature contribution, consistent with current data. 

If we introduce the dimensionlcss field <f> = p/Mp — Mp being the Planck mass — the equations of motion 
become 

(5) 



(6) 



where p and P refer to the total energy density and pressure of the standard plasma components respectively, 
and a dot denotes a derivative with respect to conformal time r\. 




III. MIMICKING RADIATION DURING BIG BANG NUCLEOSYNTHESIS 



In the standard scenario, the universe is radiation dominated during BBN. This requirement leads to theo- 
retical expectations that are in reasonable agreement with observations, at least at the 2 tr level. In particular 
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the Deuterium number density [53[ points to a value for the baryon content of the universe which is within 1-er 
equal to the independent determination of this parameter extracted from the CMB power spectrum [l[ 



Q b h 2 = 0.0223 ±0.0008 



(7) 



This result is compatible with a radiation content given by photons and three weakly interacting neutri- 
nos/antineutrinos (for a review of the present status of standard BBN see e.g. [54l. l55l. [56j ) . 

Nevertheless, present data still leave room for a non-standard energy density content, or a non standard 
value for the Hubble expansion rate. For example, the analysis carried out in [55[ shows that a conservative 
observational range for the helium-4 mass fraction constrains any extra contribution or deficit in the radiation 
energy content in the range —1.1 < AN e g < 0.8, when parameterizing it in terms of extra effective neutrino 
species. 

In view of these considerations, the main question we would like to address is whether it is still possible 
that the BD field <j> contributes during BBN in a non negligible way. The effect of </) can be easily discerned 
from Eq. ([6]). In particular, its dynamics affects the value of the gravitational constant which scales as l/</> and 
which ties the energy-density to the expansion rate. Moreover, the dynamics of <j> provides extra contributions 
to the Hubble parameter. The usual assumption is to consider a negligible role of the functions W/Wircj) (the 
"potential energy term") and K = (2lj + 3)M p (<j)/(j>) 2 /32ira 2 (the "kinetic energy term") in Eq. (JBJ) 1 , so that 
the usual BBN scenario is only possibly changed by a different value of the effective gravitational constant or, 
equivalently, of a re-scaled energy density, see e.g. [57| . The same bound previously quoted in term of N e g then 
translates into 0.9 < < 1.2. 

We now relax the hypothesis that the cf> field dynamics contributes negligibly at the epoch of BBN. Instead, 
we consider the possibility that the behavior of (f> during the temperature interval 1 McV >T> 0.01 MeV may 
be highly nontrivial and yet may still result in primordial abundances for 2 H, 3 He, 4 He and 7 Li of the correct 
order of magnitude. In particular, we consider a specific, yet interesting case, in which the universe during this 
stage is still radiation dominated, but by this we now only mean that the specific dependence of the scale factor 
a(t) oc i 1 / 2 , or equivalently in terms of the conformal time 



0(7?) = KT] 
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where k is a constant to be evaluated by cosmological observables. The condition ||5j) in turn imposes severe 
constraints on Eqs. ([5]) and ©, which under this ansatz can be recast into 



2 

-( 

'/ 



k 2 Mp (2w + 3) 
2lu + 3 



H 4tt ft 8tt 



12 



8TTG 2 



3 k 2 Mp , 



W 

167T 



1/2 



where a prime denotes a derivative with respect to a. Squaring Eq. (|10[) one obtains 

K 2 M 2 
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Substituting this expression for W into Eq. @ yields the following nonlinear ordinary differential equation for 
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Note that by writing the differential equation in this form we have excluded the case in which the field satisfies 
the first order equation 



(13) 



1 Notice that this definition of the kinetic and potential energy density for <j> is somehow imprecise, since their sum does not corre- 
spond in general to a covariantly conserved energy density. Nevertheless they are useful quantities to quantify the contribution 
of the BD field dynamics to the expansion rate and moreover, their sum is covariantly conserved in the limit of very large ui and 
negligible energy density p of standard matter which we discuss in this section, see e.g. [2Stl ■ 
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which from Eq. ([TO]) implies that the associated solution is only possible for a negative potential W. Therefore, 
we shall not consider this case further. 

Using the covariant conservation of the energy for standard matter ap' = — 3(p + P) , Eq. (JT3J) can be further 
simplified into 

„ 40 2 + 2a00' - a 2 uj(/> /2 24tt a(p + P) 



a 2 k 2 Mf> (2w + 3) 



2 T7 ■ (14) 



We are interested in a situation in which the expansion rate of the universe is of the order of the standard 
value leading to a successful BBN, and yet the contribution is relevant or even dominant. By inspection of ^ 
we see that in order for these conditions to hold, the energy density p should be suppressed by a value of > 1- 
Furthermore, the new terms, which arc dependent on the dynamics of 0, should mimic a radiation dominated 
behavior and compensate for the reduced contribution of ordinary radiation to the expansion rate. In order to 
see whether this scenario is in fact possible at all, we start by considering the simplest case, in which ordinary 
radiation and matter contribute negligibly to the kinematics. Of course, this means that we have to check a 
posteriori that, during BBN, it is possible to choose values of large enough so that p is sufficiently suppressed. 
The leading contribution to ^ will then be provided by the dynamics of 0. 

If we thus neglect the second term at the r.h.s of Eq. (fT4|) , the replacement s = <f> 1+ul linearizes the equation, 
giving 

s»= 2 -s' + A -^s , (15) 
a a z 

from which the general solution for 0(a) can be expressed as 

4 > (a)=$o\(l-C)a( 3 - A V 2 + Ca( 3+A V 2 ] 1/il+ " ) , (16) 



with A = v25 + 16 oj, $0 a normalization constant giving the value of the field at a = 1, and C determines the 
initial value of 0' via 

t^V^iJ^—.lMic-V + s] . (17) 

2(1 + uj) 

It should be noticed that in order that the solution remain valid in a neighborhood of the BBN epoch we should 
have < C < 1. Furthermore, it can be checked that the potential W fails to be positive definite if C > 1. 
There are two possible monotonic behaviors for in the two particular cases C = and C = 1 

± ( a ) = $ o a e ± , (18) 

where e± = (3 ± A)/2(l + uj) = ±8/(A =p 3). The e + and e_ indexes correspond to the field rolling towards 
larger or smaller values respectively, with an extremely slow power-law behavior since, for large uj, we have 
|e±| ~ 2/y/ui < 0.01. For all other choices of C the solution satisfies a very finely tuned condition corresponding 
to the field bouncing back at <E>0; as shown in Figure [TJ In the following therefore, we will only consider the 
two behaviors described in Eq. (fl~8"]) . Notice that since u> is very large, in the interval 0.01 < a < 100 the two 
solutions can be expressed with an accuracy better than 1% as 

0±(a) = $0 (l±^=\ogo)j . (19) 

Plugging Eq. (fTS)) into Eq. (|11[) (and neglecting p) one has 

^ = ^A4(6 + 6 £± -.4)(^)" 4A± -2^M|(i-) ±2 ^ , (20) 

where the last expression holds for large uj. 

Let us summarize the results obtained so far. By considering a potential of the form shown in Eq. (|20p . the 
expansion rate of the Universe due to a BD field behaves the same way as during a radiation dominated regime. 
Depending on the sign of the power-index in W, there are two possible behaviors for the field 0, increasing or 
decreasing with cosmic expansion. The field rolls very smoothly along the potential, and with comparable 
"kinetic" and "potential" energy densities K and W/(j> respectively. Since we derived these results in the limit 
W 3> p, we must check that this condition is compatible with the solution found. Indeed, it is easily seen that 
the term W/<fi in Eq. (|20]) is independent from the value of while the term p/0 can be made arbitrarily small 
by choosing $0 arbitrarily large. This implies that provided we choose the initial value for to be sufficiently 
large, the contribution of ordinary radiation can be neglected, and the expansion of the Universe follows a oc 7/ 
despite the fact that it is driven by a BD scalar field. 
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FIG. 1: The evolution of the 0/3>o field versus a for C = 0.5 and u = 40000. The field bounces back at a — 1 where 

4> = $0. 

IV. BBN IN A BD-DOMINATED COSMOLOGY 

In the previous section we have found that there are two solutions for which a BD field dominates cosmic 
dynamics, while keeping a radiation-dominated expansion as in Eq. (|8|). We turn now to the phenomenological 
requirement of preserving in this scenario the general agreement between BBN predictions and the observed 
light clement yields. Of course, the scenario outlined in the previous section of a pure ^-dominated universe 
cannot describe the late stages of cosmic evolution. Indeed this would be at variance with the fact that today 
we require <f> — 1. Moreover, with the potential considered in the previous section, the universe would continue 
to expand as in a radiation dominated regime until very small redshifts. We do not address these issues further 
here, but it is clear that a more complicated dynamics after the BBN period is required in order to bring 
such kind of models in accordance with "late" cosmological and astrophysical observables. In the following, in 
particular wc will assume that CMB observations arc explained as in the standard scenario. 

The set of equations governing BBN, besides equations j5|) and ([6]), is the following (see e.g. [55j ]) 

-m (21) 

-3H{p + P) (22) 
TiiXj) (23) 

3 

klLL-]. (25) 

The first two equations state respectively the conservation of the total baryon number and energy density in a 
comoving volume. The third equation is the Boltzmann equation which describes the density of each nuclide 
species, with the Ti being the rates of interaction averaged over the kinetic equilibrium distribution functions. 
The fourth equation accounts for the electric charge neutrality of the universe in terms of the electron chemical 
potential <I> e . Finally, the last equation is the Boltzmann equation for the neutrino species. 

It is clear that the field <f> never explicitly enters these equations but does appear in the evolution of the Hubble 
parameter. Therefore, as long as the evolution of H in the BD case is indistinguishable from the standard case, 
adopting initial conditions (at whatever T > 1 MeV) for the matter and radiation fields analogous to the standard 
case, the production of the light elements will also be unchanged. However, even in the scenario described in 
the previous section there is one relevant departure, since the Hubble factor evolves strictly as a -2 , while in 
the standard case this evolution is modified because e annihilation changes the number of effective relativistic 
degrees of freedom. Note that this phenomenon happens just between the two crucial temperatures for BBN; 
T ~ 0.7 MeV at which the neutron/proton ratio freezes, and T ~ 0.07 MeV when the deuterium bottleneck 
opens and light element production begins. If the BD Hubble rate matches one of the two conditions, it might 
fail to match the second one. To test if such a difference can be accommodated at all, we have run the BBN 
code described in [55l | by assuming that the dominant contribution to H is provided by the BD field, and its 
behavior is as for pure radiation as discussed in the previous section. To this end we replace the standard 
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matter/radiation energy density p in the Hubble law (and only there!) with the following paramctrization in 
terms of a single neutrino energy density 



H 



I 8tt 
3A/J, 



(26) 



where p v = 7tt 2 T i 4 /120 and T v cx a . In the limit in which the solution previously found holds, this useful 
paramctrization is exact. We ask if there is any value at all for the factor y which can accommodate at least 
4 He and 2 H abundances. Since we assume that the standard cosmology holds for a 3> obbn, we fix the baryon 
abundance to that deduced from the CMB. Eq. ([7|). Note that, since there are two conditions to fulfill with 
only one parameter, the constraint is non-trivial. We find that, for 



5-5 < y < 7.4 



(27) 



the conservative constraint on the 4 He mass fraction 0.232 
predicting a deuterium fraction 1.99 • 10~ 5 < 2 H/H < 2.66 



< Y t 
10" ; 



< 0.258 quoted in [58| is satisfied, while 
comparable with the observed one 2 /H= 



(2.78^3g)-10- 5 [53| within 2 a. For y in the range of Eq. (27]) we also found 0.97 ■ 10" 5 < 3 Hc/H < 1.07- 10" 10 , 
which is consistent with present bounds [E^, and 4.2 • 10 -10 < 7 Li/H < 5.1 ■ 10 -10 , which remains a factor ~ 3 
larger than the observed one [13, [HH . 

In terms of d>, this constraint translates into 
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or equivalently 



5.16 (T u ay < kM p Hr 



l + e± + ^ < 5.97 (T v a) 2 



(29) 



The existing bound on uj implies |e±| < 0.01, which in turn simplifies the previous constraint to the following 
one (with a percent accuracy) 



5.16 {T v a) 2 <kM p < 5.97 (T„ a) 2 



(30) 



This is the constraint on the potential W — or, more precisely, on W/<p — of Eq. ([20|) which accommodates our 
model with successful nucleosynthesis. From current cosmological measurements one would deduce T v a ~ 2 K 
and constant from the BBN epoch till now while. However, lacking a detailed modelling of the transition from 
a pure (^-dominated phase to a radiation dominated one, we do not elaborate on this bound further 2 . 

Since we have shown that BBN is compatible not only with a "minimal" effect, but also with a major role 
for a BD field, it is reasonable to expect that intermediate cases in which the contribution of cf> and of ordinary 
radiation are comparable might be viable. Of course, such scenarios generally involve a high degree of fine 
tuning in the parameters of the field 4> and/or in the choice of the potential W . Thus, for illustrative purposes, 
in the following we shall take a phenomenological approach and discuss the situation where Eq. |0 is replaced 

by 
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and again consider the case for which the BD dynamics corresponds to a pure radiation behavior, but relaxing 
the hypothesis of a vanishing contribution of ordinary radiation to the Hubble parameter. In the first equality 



2 In the "standard cosmology" , the value of k would be 

K = H a 2 = Ho afi s/Q, r , 

where ao and f2 r are the scale factor and the radiation fraction today including also neutrinos, independently of the fact that at 
least some of them are non-relativistic today. However, unless one specifies the mechanism for the field <j> to exit the epoch of 
(^-dominance, one has no way to relate the k needed at the BBN epoch with present-day cosmological quantities. 
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FIG. 2: The 1 (black), 2 (grey) and 3cr (light grey) contours in the <j> — y plane. 



of Eq. (j3Tj) we have neglected the term 4>/2acf>, which is anyway suppressed by a factor < 0.01. Running 

the modified code for this case leads to the bounds shown in Figure^ where we plot the likelihood contours 
in the (j) — y plane at 1, 2 and 3-cr level using both 2 H and 4 He experimental results as before. Two limits are 
easily recognized. For ^ > 1 we recover the case considered previously of a dominant BD contribution to H, 
and correspondingly the bound on y reported in Eq. (|27p . while for !/ « 1, we obtain the "trivial" bound on 
the effective Newton constant during BBN, namely 0.9 < <?!> < 1.2. This result shows that not only both regimes 
are viable, but also that intermediate situations are phenomenologically allowed, although the details of these 
cases are more model dependent, and the parameters in the BD sector require a high degree of fine tuning. 
It should be noticed that to impose the condition that BD behaves as radiation, or that 

W 

K + T7—L=yP»> ( 32 ) 

1O7T0 

corresponds to the choice of a particular potential. In principle one should show that this choice is possible. In 
order to see if this is indeed sensible, we use (|32p in equations ((5])-((6|) and find that a solution can be found by 
taking the time derivative of these two equations. The outcome is a system of nonlinear 2 nd order differential 
equations in a and <fi. There is only one independent initial condition — the choice of <pi — since di can be found 
by using equation and since equation §5§ becomes a cubic equation in (f>. This means that the system 
becomes uniquely determined and the solution does indeed exist and it is unique. 



V. COMMENTS AND CONCLUSIONS 



We have revisited the constraints on a Brans-Dickc scalar field, arising from the requirement of successful 
nucleosynthesis. This is an important question in light of the different roles that a Brans-Dicke field <p may play 
in various sectors of cosmology, such as inflation, baryogencsis and dark energy. 

Indeed, we have found that the BD field, i.e. Newton's constant, can have significantly different values from 
today's one, and at the same time dominate the energy density making the universe behave as if radiation dom- 
inates. The matter content is not changed at all, so that all physical quantities such as freeze-out temperatures 
and decay/production rates for elements do not differ from standard ones. 

In other words, we have found that it is possible for the BD field to mimic the effect of radiation on the 
evolution of the cosmic scale factor. In such a case is the dominant contributor to the Hubble parameter 
during BBN. We have put constraints on the BD potential using the observed light elements yields and the 
BBN predictions in such a modified scenario. This scenario requires a large initial value for > 10 so that 
the standard contribution of matter to the expansion rate is suppressed by the reduced value of the effective 
Newton constant, as well as potential for the BD field W that scales as a power law in <f>, 

W/4> cx 4> ±2 ^ . (33) 

Here oj > 40000 if we assume that the local result of the Cassini mission also holds on cosmological scales. 
We do not offer any suggestion how to explain such an high exponent. By numerically solving the set of BBN 
equations we have found that it is indeed possible to tune the overall scale of the potential W, see Eq. (f29|) . so 
as to obtain final yields for both 2 H and 4 He compatible with the observed data. 
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Of course, one may worry that large values of <p at the BBN epoch are not compatible with today's value. In 
fact, if we considered the same effective potential valid until today, the field </> would keep on being the dominant 
component of the universe, leading to a completely different evolution for the universe, e.g. an extrapolation 
would lead to value for the field (j) incompatible with today's value for the Newton gravitational constant. 
Therefore, after BBN, the potential should be chosen such that the field decreases to acceptable values, and 
the universe follows a standard evolution. The fact that today's value of (f> imposes an upper bound on the 
possibility of having very large values for 4>bbn is tantamount to saying that a fine-tuning of the choice of the 
potential is necessary in order to allow a huge change in the value of the field from BBN up to now. 

Finally, we have considered the more general scenario where both standard radiation and the BD field con- 
tribute to the expansion rate, in the case that the BD dynamics again corresponds to a pure radiation behavior. 
In this case too, by comparing the theoretical prediction with light nuclei abundances, we have found that, 
depending on the value of the BD field during BBN and its contribution to H, there is a region in this two- 
parameter space leading to successful nucleosynthesis. This region shown in Figure [2] extends from the high </> 
regime already discussed down to the standard scenario with 0=1 showing that provided the <fi potential is 
suitably tuned, basically all values of the BD field can be shown to be compatible with the abundances of both 
2 H and 4 He produced during BBN. 

An interesting feature of the scenario considered in this paper is that one is basically free to adjust the 
time-temperature relation, which analytically means rescaling the initial conditions for the thermodynamical 
quantities of the fluids in the plasma. We have limited our discussion to the case where H (a) around the BBN 
epoch has the same evolution it would have in a radiation-dominated plasma, keeping the particle densities fixed 
at the same value of the standard case as our initial conditions. This is not a fine-tuning required a priori, but 
an useful working hypothesis which authomatically ensures that the rates for the scattering processes (including 
the nuclear reactions) are unchanged. However, one may relax this assumption, and thus the allowed region in 
the model space may enlarge considerably, to the price of major (and model-dependent) modifications in the 
physical processes in the plasma. We do not discuss further this point in our phenomenological approach, since 
our purpose in this article was to show that exotic scenarios where the BD field is driving the expansion rate 
at the BBN epoch are indeed possible. Nonetheless, this feature should be kept in mind when characterizing 
specific models (as e.g. quintessence potentials with solutions tracking radiation) with a cosmic evolution before 
and after the BBN epoch determined self-consistently within the theory. 
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